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Abstract 

Classical approach of solvability problem has shed much light on what we can solve and what 
we cannot solve mathematically. Starting with quadratic equation, we know that we can solve 
it by the quadratic formula which uses square root. Polynomial is a generalization of quadratic 
equation. If we define solvability by using only square roots, cube roots etc, then polynomials 
are not solvable by radicals (square root, cube root etc). We can classify polynomials into simple 
(solvable by radicals) and complex (not solvable by radicals). We will use the same metaphor 
to separate what is solvable (simple part) and what is not solvable (complex part). 

This paper is a result of our presentation at a University of Houston seminar. In this paper, 
we will study seismic complexity through the eyes of solvability. We will investigate model 
complexity, data complexity and operator complexity. Model complexity is demonstrated by 
multiple scattering in a complex model like Cantor layers. Data complexity is studied through 
Betti numbers (topology /cohomology). Data can be decomposed as simple part and complex 
part. The simple part is solvable as an inverse problem. The complex part could be studied 
qualitatively by topological method like Betti numbers. Operator complexity is viewed through 
semigroup theory, specifically through idempotents (opposite of group theory). Operators that 
form a group are invertible (solvable) while semigroup of operators is not invertible (not solvable) 
in general. 
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Short history of solvability: 
« simplicity » 

Quadratic equation (Complex numbers) 
Fundamental theorem of algebra (Gauss) 
Solvable by radicals (Abel) 

Group theory (Galois/solvable group, Lie/continuous group) 
Operator theory (Hilbert, infinite dimensional space) 

« complexity » 

Numerical complexity (Poincare /chaos) 

Model complexity (Cantor / layers) 

Data complexity (Betti number/ diffractions) 

Operator complexity (Idempotent semigroup/ 2-way wave ) 



Complex Decomposition (Data or operator) 

= Simple part (quantitative/solvable) 

+ Complex part (qualitative/interpretable/not noise) 
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Presentation based primarily on 4 papers published on arXiv 
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Part 1 : Model complexity (Cantor / layers): 
Cantor set models and solvability of multiples 
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The Cantor Set 

Construction: repeatedly removing middle thirds of line segments. 

The Cantor set = all points that are not removed at any step in this infinite 
process. 



rj The Cantor set, in six iterations i 
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Example 1 - Discrete Horizontal Cantor Set Model: C4 



Approximate 1 st -order 

Model (Vp=Constant) P Multiples (M) = P * P 




Figure 6: Slide 6 



4 



C is sparse and totally disconnected 
C+C=l is dense and nowhere disconnected 

Multiples may not be solvable for a Cantor set medium 

The Cantor Set Models: No Means + No Variance 
(No Effective Medium) 
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Example 2 - Discrete Vertical Cantor Set Model: VC6 
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Example 2 - Discrete Vertical Cantor Set Model: VC6 
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Near and Far Stacks See Different Models of Complex Topology 
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Example 3 - Sierpinski Carpet: 
a 2-D Generalization of the Cantor Set 



Construction: divide a square into 9 equal sub-squares in a 3-by-3 grid, and 
remove the center sub-square; repeat the process for the remaining 8 sub- 
squares; and apply the process recursively. 
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Example 3 - Discrete Sierpinski Carpet Model: SP4 
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Example 3 - Discrete Sierpinski Carpet Model: SP4 




Figure 13: Slide 13 

Example 3 - Complex topology in the overburden 

Model (V=Constant) Modeled CMP gathers 
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Results are generated from a deterministic model without adding "random 
noise" 
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Example 3 - Sierpinski Carpet Models: SP4 



— CMP-sorted Gathers 




Flat reflector 
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Model complexity Summary 

The Cantor set: 

- The Cantor set C is sparse and totally disconnected; while C+C=l is 
a continuum that is dense and nowhere disconnected 

- The Sierpinski Carpet is a 2-D generalization of the Cantor set. 
Example 1 (horizontal Cantor set models): 

- Multiples are not solvable for the Cantor set model 

- No means + no variance (no effective medium) 
Example 2 (vertical Cantor set models): 

- "Noise" from inversion might be an indication on the earth's complex 
topology 

Example 3 (Sierpinski Carpet models): 

- A model with complex topology can generate realistically looking, 
complex wave field 
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Part 2: Data complexity (Betti number/ 
diffractions): 
Geometric simplicity as a migration criterion 
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Betti numbers and homology groups 

• Homology is a tool to measure connectivity and holes in a 
topological space; Betti numbers are dimensions of 
homological groups. 

• The nth Betti number, Bn, represent "n-dimensional holes" 

• For a 2D image, only BO and B1 are non-trivial, Bn=0 (n>1) 

• The th Betti number, BO, counts the number of connected 
components 

• The 1 st Betti number, B1 , counts the number of "holes" 

• The following 3 images are equivalent homologically, i.e. 1 
connected component and 1 hole in each 




B0=1 B1=1 
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BO & B1 : more examples 



• 


• • 


Figure 1a. B0=1 B1=0 


Figure 1b. B0=2 B1=0 
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Figure 1c. B0=1 B1=1 


Figure 1d B0=1 B1=3 
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Betti numbers and the 
Computational Homology Project (CHOMP) 



• CHOMP is a US-based, international open- 
source development project for promoting the 
application of computational topology to global 
analysis of nonlinear spaces and functions. 

• Web address 

• http://chomp.rutgers.edu/ 

• We take advantage of the latest developments in 
CHOMP to carry out efficient computation of 
Betti numbers on seismic data 
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Betti numbers, geometric complexity and wave propagation 



Increasing geometric complexity (B1 grows) 




Wave forward propagation {modeling or de-migration) 



Wave backward propagation (focusing or migration) 
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Figure 2: A synthetic dataset containing three diffr actors. 
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Figure 3: Zoom-in view on the semigroup migration. 
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Migration velocity 



Figure 4: Betti number Bl as a function of migration velocities for the syn- 
thetic dataset in Figures 2 and 3. 
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Figure 5: A 2D stack section before migration. 
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Figure 6: Zoom-in view on a portion of the dataset from the semigroup 
migration. 
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Figure 7: Betti number Bl versus migration velocity, corresponding to Fig 
6. that reaches minumum at a migration velocity of around 1800 m/s. 
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Examples of migration-operator 
related "noise" 




Primaries interfered by numerical artifacts 
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Proposed steps for reducing migration noise 

• Forward model or de-migrate 

• Apply desired process, e.g., mixing, to 
attenuate noise, with one parameter 

• Compute Betti numbers on time slices to 
infer level of geometric complexity 

• Analyze Betti numbers and determine 
optimal processing parameter 

• Inverse model or re-migrate 
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Example 1 : scatterers 




Figure 2. Input scatterers form a disordered system 
(it is not coherent and not random) 
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Example 1 : forward model or de-migration 
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Example 1 : migration with 3% error 
we desire to re-focus to minimize the migration "smiles" 



r 




Figure 4. Migration output with 3% velocity error 
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Example 1 : compute Betti numbers for a varying 
processing parameter (mixing radius in this example) 



Figure 5a. Diffraction pattern after 7x7 mix 



Figure 5b. Timeslice of diffraction pattern after 7x7 mix 



Figure 33: Slide 33 



Example 1 : reducing geometric complexity defined by B1 



Mix 7x7 shows rapid decline in 
Betti number B1 (green curve) 




a 5 IB 15 20 

Hixing radius <# points) 

Figure 6. X-axis is amount of mixing and y-axis is Betti numbers 



Complex 



Less complex 



Figure 3b. Timeslice of diffraction pattern 



Figure 5b. Timeslice of diffraction pattern after 7x7 mix 
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Example 1 : re-migrated with 
geometric simplicity 




Figure 4. Migration output with 3% velocity error Figure 7. Migrated output after 7x7 mi: 
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Example 2: real data 




Primaries interfered by numerical artifacts 
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Example 2: re-migrated with 
geometric simplicity 



Original Reduced complexity 
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Geometric simplicity and seismic imaging 

• Geometric simplicity as defined by the Betti 
numbers can be used for quantitatively 
determining impact of processing parameters 

• Seismic imaging should always start from simple 
geometry of the input data to build intuition for 
the migration operator 

• As we build better understanding of the data and 
the operator, we can increase the Betti numbers 
and increase the geometric complexity for 
migration 
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Summary 

• Geometric simplicity can be defined by 
Betti numbers 

• CHOMP made it possible for efficient 
computation of Betti numbers 

• Examples demonstrated usage of Betti 
numbers for reducing migration artifacts 

• This approach can be easily generalized 
to other applications, such as tomography 
and waveform inversion 
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Part 3: Operator complexity (idempotent 
semigroup) : 

2-way wave field extrapolation and imaging 
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Complex decomposition of operators: 

An example of wave field extrapolation of 1 D acoustic wave equation 
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Diagonalize: Another way of further analyzing Eq,(4) is to examine eigenvalues and eigenvectors of A, i.e., 
A = VBV -1 . where V is the eigenvector matrix. D diagonal matrix of the eigenvalues. 
Eq.(4) becomes 



where the diagonal eigenvalue matrix is 

D = 



-i- 



Apply projection: g = V f 



(6) 



(7) 



Op.ra.or: L = ( = Pi + P 2 = ei<^+ e 2 < ^ 

where 
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and C2 = 
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are idenipotents which correspond to upward-propagating and 
downward-propagating components of the wavefield, respectively 
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Complex decomposition of operators: Methodology 



We will first introduce some mathematical concepts in semigroup and group theory. In order to 
make it less abstract, we will think of operators as matrices since matrices are a good surrogate 
to discuss abstract ideas. If we take all n X n matrices under matrix multiplication as a whole 
system, it is a semigroup and it has the semigroup property of associative law A(BC) = (AB)C. 
This system includes both invertible and noninvertible matrices. We will come back to the idea of 
semigroup and idempotents in semigroups. 

If we restrict a subset to all invertible matrices, then it is a group since each n x n invertible 
matrix A by definition will have a B so that AB = BA — I the identity matrix. Two matrices are 
equivalent if A = UBV where U and V are invertible. A desirable property is to write a matrix as 
A = UDV where D is a. diagonal matrix. More generally, it is desirable to have A = UPV where 
P is an idempotent with PP = P. 
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Multiplication Zero Addition 

(less energy) (more energy) 



Figure 1: A lower energy, stable state can be reached by multiplication of idempotents. 
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Complex decomposition of operators 



Decomposition of operators (with predetermined basis) 

The simplest way to decompose an operator is to use a predetermined basis like Fourier basis. The 
advantage of predetermined basis is to convert a differential equation into an algebraic equation 
with orthogonal basis like fk (/ is frequency and k is wavenumber). The evanescent wave is just 
filtering certain dips (i.e. certain fk numbers) (Figure 2). The disadvantage is the difficulty of 
handling large velocity variations. 

Any predetermined basis forces an operator decomposition into predetermined geometry like 
Fourier basis. The predetermined basis maps the operators into fixed shapes given by the basis 
elements. 

Decomposition of operators (without predetermined basis) 

The second method is to lump the derivatives with the velocity field so it is more sensitive to the 
velocity variations. Finite difference and Kirchhoff extrapolation are such methods but they can 
have practical limitations for frequencies and dips. Integral methods like Kirchhoff extrapolation 
(see Berryhill) are more geometrically motivated. An excellent numerical treatment can be found 
in Sandberg and Beylkin 2009 where positive eigenvalues are viewed as evanescent energy. The 
negative and zero eigenvalues are kept as propagating waves in the extrapolation operator. This is 
another way to separate the wavefield into simple part and complex part. 
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Complex decomposition of operators 

Conclusion 

Solvability by semigroup is a way to use noninvertible operators to approximate the original oper- 
ator. We have chosen in this paper to use idempotents P where PP = P. Semigroup theory is not 
restricted to idempotents but idempotents are easier to demonstrate since they resemble diagonal 
matrices (see Davies. Wilansky). They are also familiar to physicists in terms of density matrix 
with or 1 eigenvalues (see Bowler et al) viewed as basic building blocks. 

Multiplication of noninvertible operators becomes less energetic or contractive or lossy or more 
limited. Addition of noninvertible operators becomes more energetic or expansive or generating 
larger operator. Toggling between multiplication and addition of the idempotents gives us the 
stable part (multiplicative) and the generative part (addition). This is the essence of downward 
continuation in wavefield extrapolation. The schematic is captured with the diagram in the end 
(Figure 1). 

The simple part of the wavefield extrapolation operator is written as the sum of idempotent 
operators ( A = PI + P2 4- P3 + . . . ) where A is the operator and Pj is an idempotent with 
PjPj — Pj and PjPk — if Pj and are different. We dropped the phrase equivalent to idempotent 
to emphasize the importance of idempotents. The simple part of the operator is the propagating 
operator to extrapolate the wavefield. It can be studied through semigroup of operators. The 
complex part of the operator is the evanescent energy which should not be propagated. The 
complex part needs qualitative mathematics like topology to classify and understand. 
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Short history of solvability: 
« simplicity » 

Quadratic equation (Complex numbers) 
Fundamental theorem of algebra (Gauss) 
Solvable by radicals (Abel) 

Group theory (Galois/solvable group, Lie/continuous group) 
Operator theory (Hilbert, infinite dimensional space) 

« complexity » 

Numerical complexity (Poincare /chaos) 

Model complexity (Cantor / layers) 

Data complexity (Betti number/ diffractions) 

Operator complexity (Idempotent semigroup/ 2-way wave ) 

Complex Decomposition (data is operator and operator is data) 
= Simple part (solvable) + Complex part (interpretable) 

Toy Modeling : Simple + Simple = Complex (C+C=l with Cantor + Cantor = 

Continuum) 

Complex Decomp of data : Complex = Simple + Complex (low Betti + high Betti) 
Complex Decomp of operator : Simple + Simple = Complex (low rank + low rank = high 
rank) 
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